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STUDIES  IN  SATURATED  POOL  BOILING  OF  WATER 
Technical  Note  N-1209 
ZFXX.512.001.004 
by 

S.  C.  Garg 


ABSTRACT 

A  photographic  study  of  saturated  nucleate  pool  boiling  of  water 
on  horizontally  mounted  wires  and  tubes  was  carried  out.  The  effects 
of  the  diameter  of  the  heating  surface,  system  pressure  and  heat  flux 
upon  the  bubble  diameter  at  departure  from  the  heating  surface  were 
investigated.  Five  platinum  wires  of  diameter  between  0.001  and 
0.026  inch  and  two  platinum  tubes  of  diameter  0.05  and  0.10  inch  were 
used  at  pressures  between  1.0  and  14.7  psia.  Departure  diameter  of 
more  than  1,800  discrete  bubbles  was  measured  at  heat  fluxes  up  to 
125,000  Btu/hr-ft^  by  a  frame  by  frame  analysis  of  6,100  feet  of 
16  mm  film  taken  at  3,000  and  5,000  frames  per  second. 

The  bubble  diameter  data  were  compared  with  analytical  expres¬ 
sions  by  Owens,  Ivey  and  Morris  and  Garg.  An  empirical  expression 
is  proposed  which  successfully  correlates  the  bubble  diameters  in 
saturated  nucleate  pool  boiling  of  water  at  all  pressures.  Recommen¬ 
dations  are  made  to  systematically  investigate  the  effects  of  surface 
roughness  and  the  degree  of  subcooling  upon  the  bubble  diameters  at 
departure  in  nucleate  pool^jjdiTing . 
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A  photographic  study  of  saturated  nucleate  pool  boiling  of  water 
on  horizontally  mounted  wires  and  tubes  was  carried  out.  The  effects 
of  the  diameter  of  the  heating  surface,  system  pressure  and  heat  flux 
upon  the  bubble  diameter  at  departure  from  the  heating  surface  were 
investigated.  Five  platinum  wires  of  diameter  between  0.001  and 
0.026  inch  and  two  platinum  tubes  of.  diameter  0.05  and  0.10  inch  were 
used  at  pressures  between  1.0  and  14.7  psia.  Departure  diameter  of 
more  than  1,800  discrete  bubbles  was  measured  at  heat  fluxes  up  to 
125,000  Btu/hr-ft^  by  a  frame  by  frame  analysis  of  6,100  feet  of 
16  mm  film  taken  at  3,000  and  5,000  frames  per  second. 

The  bubble  diameter  data  were  compared  with  analytical  expres¬ 
sions  by  Owens,  Ivey  and  Morris  and  Garg.  An  empirical  expression 
is  proposed  which  successfully  correlates  the  bubble  diameters  in 
saturated  nucleate  pool  boiling  of  water  at  all  pressures.  Recommen¬ 
dations  are  made  to  systematically  investigate  the  effects  of  surface 
roughness  and  the  degree  of  subcooling  upon  the  bubble  diameters  at 
departure  in  nucleate  pool  boiling. 
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INTRODUCTION 


Boiling  of  a  liquid  on  a  metal  surface  has  long  been  recognized 
as  one  of  the  best  methods  of  obtaining  very  high  heat  transfer  rates 
for  only  a  few  degrees  temperature  difference  between  the  heating  sur¬ 
face  and  the  liquid.  Consequently,  this  method  of  heat  transfer  pre¬ 
dominates  in  most  high  power  density  systems.  Boilers,  desalination 
equipments,  heat  exchangers,  boiling  water  nuclear  reactors,  and  cooling 
of  high  power  density  electronic  equipments  are  some  of  the  better 
known  applications  of  the  boiling  phenomenon.  In  systems  like  these 
and  others,  it  is  desirable  to  obtain  as  high  a  heat  transfer  coefficient 
as  possible  for  optimum  performance  and  for  economic  reasons.  A 
limitation  to  the  power  density  is,  however,  reached  when  there  is  a 
sufficient  bubble  density  on  the  heating  surface  to  cause  a  transition 
to  the  partial  or  total  film  boiling  region. 

The  importance  of  bubble  formation,  growth,  and  departure  size  in 
determining  the  heat  transfer  characteristics  of  a  heating  surface  has 
been  discussed  in  detail  in  a  recent  literature  survey  by  Garg.^  In 
heat  pipes  employing  mesh  screen  wicks,  if  the  spacing  between  con¬ 
secutive,  small  diameter  wires  was  larger  than  the  bubble  diameter  at 
departure  under  the  heat  pipe  operating  conditions,  there  is  a  possibil¬ 
ity  that  the  vapor  may  not  be  trapped  in  the  wick  thereby  eliminating 
early  dryout  as  a  limitation  of  the  heat  pipe  performance. 

It  was,  therefore,  felt  that  an  understanding  of  the  effect  of  the 
heating  element  diameter  upon  bubble  diameter  at  departure  over  wide 
ranges  of  pressure  and  heat  flux  may  add  to  the  improved  performance  of 
heat  pipes  and  other  equipments  utilizing  boiling  heat  transfer.  Con¬ 
sequently,  an  experimental  investigation  of  saturated  nucleate  pool 
boiling  was  carried  out  tc  systematically  determine  the  effects  of 
diameter  of  the  heating  surface,  pressure  and  heat  flux  upon  the  bubble 
diameters  at  departure.  The  bubble  diameter  da  ;a  obtained  during  this 
investigation  were  compared  with  correlations  found  in  the  literature. 

A  new  empirical  correlation  is  proposed  which  correlates  the  available 
bubble  diameter  data  for  water  in  saturated  pool  boiling  in  the  discrete 
bubble  region  at  all  pressures. 


DESCRIPTION  OF  THE  EXPERIMENTAL  APPARATUS 


Test  Section 

Heating  Element.  Electrical  heating  was  considered  as  the  most 
suitable  means  of  providing  high  heat  fluxes  on  the  surface  of  small 
diameter  wires  and  tubes  for  this  investigation.  Besides  the  fact 
that  there  is  no  other  practical  way  of  heating  wires  of  diameter  as 
small  as  0.001  inch,  there  are  a  number  of  advantages  in  using  electri¬ 
cal  heating.  Among  them  are: 

a.  Heat  flux  is  uniform. 

b.  High  heat  flux  densities  can  be  easily  obtained. 

c.  Control  and  measurement  of  heat  input  rates  are  relatively 
simple  and  accurate. 

d.  The  wire  or  tube  average  outer  surface  temperature  can  be 
easily  calculated  from  the  electrical  measurements  if  the  material  of 
test  elements  has  a  high  enough  temperature  coefficient  of  electrical 
resistivity. 

Furthermore,  spot  welding  a  thermocouple  on  the  outer  surface  of  a 
small  wire  or  tube  to  measure  its  temperature  can  cause  serious  distor¬ 
tions  in  the  characteristics  of  the  local  surface  in  causing  nucleate 
boiling  and  in  its  temperature.  For  these  reasons,  it  was  decided  to 
use  heating  elements  made  of  platinum.  The  temperature  coefficient  of 
electrical  resistivity  of  platinum  of  0.002  per  degree  F  was  considered 
high  enough  to  permit  accurate  calculation  of  the  outer  wall  temperature 
of  the  heating  elements. 

Seven  different  diameters  were  selected  for  investigating  the 
effect  of  diameter  of  the  heating  surface  upon  the  bubble  diameter  at 
departure  in  saturated  nucleate  pool  boiling;  viz,  0.001,  0.002,  0.005, 
0.010,  0.025,  0.050,  and  0.100  inch.  However,  the  low  resistances  of 
0.05-  and  0.10-  inch-diameter  platinum  wires  requiring  high  input 
currents  to  obtain  the  heat  flux  levels  commensurate  with  nucleate 
boiling  investigations  necessitated  their  replacement  by  platinum  tubes 
of  equal  outer  diameters. 

Since  the  effect  of  surface  finish  upon  the  bubble  diameter  at 
departure  was  not  the  subject  of  this  investigation,  only  smoothly 
drawn  platinum  wires  and  tubes  were  used.  No  effort  was  made  to 
examine  their  surface  condition.  To  minimize  end  effects  due  to 
conduction,  it  was  decided  to  use  approximately  6-inch  heated  lengths 
of  these  elements.  The  length  to  diameter  ratio,  therefore,  varied 
between  60  and  6,000.  To  minimize  the  effect  of  disturbances  caused  at 
a  location  on  the  heating  surface  by  bubbles  departing  from  neighboring 
locations,  the  heating  elements  were  mounted  horizontally. 
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Electrical  Conductors.  The  conductors  were  fabricated  from  3/4-inch 
square  section  brass  bars  to  cause  negligible  heat  generation  at  the 
estimated  maximum  current  of  about  100  amperes,  and  chrome-plated  to 
prevent  contamination  of  water  in  the  test  tank.  The  heating  elements 
were  clamped  to  the  conductors  using  specially  designed  fittings,  shown 
in  Figure  1. 

Test  Tank.  For  success  in  high  speed  photographic  work,  it  was 
necessary  to  obtain  as  much  transmission  of  light  to  the  heating  surface 
as  possible.  It  was  therefore  decided  to  design  a  completely  transparent 
tank.  Glass  and  acrylic  were  considered.  Acrylic  was  chosen  because 
of  the  ease  with  which  it  can  be  machined.  Based  upon  strength  and 
deflection  calculations,  the  required  thickness  of  acrylic  sheets  from 
which  the  tank  was  fabricated  was  determined  to  be  3/4  inch. 

A  chrome-plated  brass  frame  of  outside  dimensions  11"  x  10"  x  6" 
was  fabricated  by  silver-brazing  1-inch-square  bars  at  the  corners 
and  all  the  outside  faces  were  machined  and  grooved  to  retain  0.125- 
inch-diameter  rubber  gaskets.  To  minimize  leaks,  machined  acrylic 
plates  were  screwed  to  the  frame  to  form  the  test  tank,  leak  proof 
seals  being  formed  by  Buna-N  rubber  gaskets  on  all  sides. 

Penetrations  for  leads  to  measure  voltage  drop  across  the  heating 
element,  a  thermocouple  to  measure  the  temperature  of  water  in  the  tank, 
brass  busbars,  and  connection  to  a  manometer  were  provided  through  the 
top  plate.  A  drain  valve  was  connected  to  the  tank  through  the  bottom 
plate.  Penetrations  for  an  auxiliary  heater  located  near  the  bottom  of 
the  tank  to  control  the  bulk  liquid  temperature  were  provided  in  the 
back  plate.  Finally,  removal  of  vapor  generated  during  boiling  and 
return  of  the  condensate  were  achieved  through  a  1^-inch-diameter  bore 
through  the  right  side  of  the  tank. 

Vacuum  System 

The  pressure  range  selected  for  this  investigation  was  from  1.0  to 
14.7  psia.  The  reasons  for  choosing  this  range  were  primarily  experi¬ 
mental,  namely, 

a.  Bubbles  are  large  at  low  pressures.  The  effect  of  heater 
diameter  upon  the  bubble  diameter  can  be  investigated  with  relative  ease. 

b.  Temperatures  are  low  which  makes  it  possible  to  use 
acrylic  for  fabrication  of  tank,  and 

c.  There  are  large  changes  in  fluid  properties  over  this 
pressure  range.  For  instance  there  is  a  2.4  fold  change  in  viscosity, 
a  1.2  fold  variation  in  surface  tension,  a  12.5  fold  variation  in  vapor 
density,  and  2.5  fold  variation  in  Prandtl  number. 


3 


Vacuum  was  created  and  maintained  by  a  vacuum  pump  connected  to 
the  tank  through  the  1^- inch-diameter  bore  connection  via  a  water- 
cooled  condenser,  which  condensed  most  of  the  vapor  generated  during 
boiling  and  returned  the  condensate  to  the  tank,  and  controlled  by 
adjusting  controlled-leak  needle-valves  located  downstream  from  the 
condenser.  A  mercury  manometer  connected  through  the  top  plate  of 
the  tank  provided  accurate  measurement  of  vacuum  in  the  vapor  space 
of  the  tank. 

Power  Supply 

To  avoid  any  influence  of  the  50  cycles  per  second  current  and 
voltage  ripples  on  bubble  initiation,  growth  or  departure,  it  was 
decided  to  provide  dc  heating  to  the  various  heating  elements.  A 
Sorensen  solid-state  dc  power  supply,  Model  DCR20-250A,  was  used  to 
provide  the  needed  electrical  power.  The  power  unit  was  equipped  witli 
controls  to  provide  continuous  variation  of  voltage  between  0  and  20 
volts  and  that  of  current  between  0  and  250  amperes. 

A  chrome-plated,  inconel-sheathed,  mineral-insulated  resistance 
element  of  about  3  kw  capacity  was  located  at  the  bottom  of  the  tank 
and  used  to  raise  the  temperature  of  water  in  the  test  tank  to  the 
desired  level  in  a  short  time. 

Photographic  Technique 

Camera.  A  Fastex  camera  capable  of  up  to  8,000  frames  per  second 
was  used  in  this  investigation  with  100-foot  rolls  of  16  mm  Eastman-4X 
high  speed  film. 

Since  bubble  nucleation  is  random  in  both  time  and  location  on  the 
heating  surface,  the  choice  of  both  the  field  of  view  and  rate  of  boil¬ 
ing  had  to  be  a  compromise.  The  field  of  view  was  largely  determined 
by  the  departure  size  of  bubbles.  The  bubble  generation  rate  in  the 
field  of  view  had  to  be  high  enough  to  ensure  capturing  some  bubbles 
during  the  short  exposure  time  of  the  100-foot  rolls,  and  low  enough 
to  prevent  gross  conglomeration  of  bubbles  and  departure  from  nucleate 
boiling  causing  burnout  of  the  heating  element.  Depending  upon  the 
pressure  in  the  tank,  an  exposure  speed  of  3,000  or  5,000  frames  per 
second  was  used. 

Lighting .  The  lighting  system  consisted  of  two  1000-watt  diffused 
lights  located  at  the  back  of  the  tank  and  two  600-watt  incandescent 
lights  located  in  fruit  of  the  tank.  All  four  lights  were  at  a  45- 
degree  angle  with  the  camera  axis  and  about  10  inches  away  from  the 
canter  of  the  front  or  back  wall  of  the  tank.  On  the  back  tankwall, 
in  the  field  of  view  of  the  camera,  a  white  ground  glass  was  mounted 
to  provide  adequate  contrast  for  filming.  The  infrared  radiation  from 
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Che  light  sources  was  filtered  by  placing  0.5-inch-thick  transparent 
tanks  containing  5%  copper  sulphate  solution  immediately  in  front  of 
the  light  sources.  The  aperture  depended  upon  the  exposure  speed, 
the  distance  between  camera  and  the  heating  element,  and  the  lens 
used.  The  arrangements  of  lights,  camera,  and  infrared  filters  in 
relation  to  the  tank  are  shown  in  Figure  2.  A  schematic  of  the  ex¬ 
perimental  arrangement  is  shorn  in  Figure  3.  Photographs  of  the 
experimental  arrangement  showing  the  various  components  are  given  in 
Figures  4  and  5. 

METHOD  OF  MEASUREMENT 

Heat  Flux.  To  permit  calculation  of  the  heat  flux  on  the  outer 
surfaces  of  the  various  heating  elements,  it  was  necessary  to  accur¬ 
ately  measure  their  inner  and  outer  diameters,  the  heated  lengths,  and 
the  voltage  drop  across  and  the  current  flowing  through  them.  Details 
of  the  techniques  used  to  measure  the  diameters  and  lengths  of  the 
various  heating  elements  are  given  in  Appendices  A  and  B.  Voltage 
drop  across  the  heated  length  of  each  element  was  measured  by  a  Leeds 
and  Northrup,  Model  K-4,  potentiometer  through  a  John  Fluke's  Model 
725A  voltage  divider.  The  voltage  drop  across  a  known  resistance  in 
series  with  the  heating  element  provided  the  value  of  current  flowing 
through  it.  Depending  upon  the  current  flowing  through  the  heating 
element,  one  of  the  four  calibrated  shunt  resistances  was  used  for 
this  measurement.  The  four  shunts  were  rated  at  25,  50,  100,  and  300 
amperes  and  their  respective  resistances  were  0.002024,  0.001000, 
0.00050056,  and  0.00026667  ohm. 

The  resistivity  of  platinum  is  strongly  dependent  upon  its  purity. 
Therefore,  it  was  felt  necessary  to  carry  out  resistivity  and  tempera¬ 
ture  coefficient  of  resistivity  measurements  for  each  of  the  seven 
heating  elements.  The  procedure  used  for  these  measurements  and  the 
results  obtained  are  detailed  in  Appendix  B. 

Temperature.  The  temperature  of  water  in  the  test  tank  was  measured 
by  a  calibrated  copper-constantan  thermocouple  located  at  the  heating 
element  level  but  at  a  distance  of  1^-inch  away  from  it  in  the  horizontal 
plane.  The  average  outer  wall  temperature  of  the  heating  element  was 
determined  indirectly  from  the  voltage  drop  across  and  the  current 
flowing  through  it.  Allowance  was  made  for  temperature  variation  across 
the  cross-section  of  the  heating  element.  The  procedure  is  described  in 
detail  in  Appendix  C. 

Pressure .  As  described  earlier,  vacuum  in  the  vapor  space  of  the 
test  tank  was  measured  by  a  mercury  manometer.  Pressure  at  the  heating 
surface  was  then  determined  from  values  of  the  barometric  pressure,  the 
height  of  water  above  the  heating  surface,  and  vacuum  in  the  vapor  space. 
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Bubble  Diameter  at  Departure.  Frame  by  frame  analysis  of  the  high 
speed  films  was  carried  out  using  a  Vanguard,  Model  C-11D,  motion 
picture  analyser.  Bubble  diameters  at  departure  from  a  heating  surface 
were  measured  for  either  all  bubbles  appearing  on  a  film  or  a  limited 
number  of  bubbles  considered  representative  of  the  film  if  there  was 
a  large  number  of  bubbles  on  the  film.  Magnification  of  the  actual 
bubble  size  on  the  projected  image  in  the  analyser  was  determined  from 
fiducial  marks  made  at  a  known  distance  apart  on  a  stainless  steel 
strip  which  was  located  in  the  field  of  view  of  the  camera  in  the 
vertical  plane  of  the  heating  element.  For  each  high  speed  film  taken, 
the  system  pressure,  the  heat  flux  on  the  outer  surface  of  the  heating 
element,  the  average  wall  to  saturation  temperature  difference,  the 
number  of  bubbles  measured,  and  the  maximum,  minimum,  and  average 
departure  diameters  of  all  bubbles  measured  were  recorded  separately 
for  each  heating  element. 

OPERATING  PROCEDURE 

Assembly  and  Cleaning  Procedure 

The  conductors,  the  platinum  test  element,  the  bulk  water  temper¬ 
ature  measuring  thermocouple,  and  leads  for  measuring  voltage  drop 
across  the  heating  element  were  appropriately  assembled  and  connected 
through  the  top  plate  of  the  tank.  The  auxiliary  heater  was  located 
at  the  bottom  of  the  test  tank  and  connected  to  the  power  supply 
through  the  rear  plate  of  the  tank.  All  the  tank  faces  except  the 
front  face  were  next  screwed  to  the  brass  frame,  leak  proof  joints 
being  obtained  by  rubber  gaskets  on  all  sides  and  in  all  penetrations 
through  the  tank  walls.  The  assembly  was  then  thoroughly  cleaned  with 
acetone  to  remove  traces  of  dirt  and  oil  followed  by  several  rinses 
with  demineralized  water.  The  front  face  was  next  screwed  to  the 
brass  frame  and  the  test  tank  connected  to  the  vacuum  system  through 
the  glass  condenser. 

The  tank  was  filled  with  demineralized  water  to  a  height  of  about 
3^  inches  above  the  heating  element,  and  the  water  heated  to  about  105°F 
using  the  auxiliary  heater.  Pressure  in  the  tank  was  then  reduced  until 
the  temperature  of  water  corresponded  to  the  saturation  temperature  at 
the  pressure  in  the  tank.  Boiling  was  then  started  on  the  auxiliary 
heater  and  maintained  for  a  few  hours.  The  vapor  generated  during 
boiling  was  condensed  in  the  condenser  and  returned  to  the  test  tank. 
This  cleaning  procedure  was  carried  out  to  remove  traces  of  dirt  from 
the  glass  condenser  which  might  enter  the  tank  during  testing.  After 
several  hours  of  boiling,  the  vacuum  in  the  tank  was  broken  and  the 
water  drained.  The  tank  was  next  thoroughly  rinsed  and  filled  with 
demineralized  water.  This  cleaning  procedure  was  repeated  several 
times  to  ensure  a  thoroughly  clean  system  at  the  commencement  of  tests. 
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Pre-test  Procedure 

The  test  tank  was  filled  with  demineralized  water  to  a  height  of 
about  four  inches  above  the  heating  element,  and  heated  by  auxiliary 
heater  to  a  temperature  some  ^ -5°F  above  the  saturation  temperature 
corresponding  to  half  the  test  pressure.  The  pressure  in  the  test  tank 
was  then  reduced  gradually  to  about  half  the  test  pressure,  thereby 
causing  bulk  boiling.  Water  was  then  boiled  on  the  heating  element 
for  about  half  an  hour  under  saturated  conditions  to  minimize  the 
quantity  of  gas  in  the  bulk  liquid  and  that  absorbed  on  the  heating 
surface.  Vacuum  in  the  test  tank  was  then  broken  and  the  water 
temperature  raised  to  the  saturation  temperature  for  the  test  pressure 
by  increasing  the  voltage  drop  across  the  auxiliary  heater. 

Test  Procedure 

Pressure  in  the  vapor  space  of  the  test  tank  was  reduced  slowly 
until  pressure  at  the  heating  surface  corresponded  to  the  test  pressure. 

A  small  dc  voltage  was  applied  across  the  test  element  and  increased 
slowly  until  boiling  was  observed  on  its  surface.  Power  input  to  the 
element  was  then  adjusted  to  obtain  the  desired  rate  of  bubble  forma¬ 
tion.  When  steady  state  conditions  were  reached,  values  of  current 
flowing  through  and  the  voltage  drop  across  the  heating  element,  vacuum 
in  the  vapor  space  above  the  heating  element,  barometric  pressure, 
height  of  water  above  the  heating  surface,  and  bulk  temperature  of 
water  were  measured. 

For  high  speed  photography,  the  power  to  auxiliary  heater  was 
turned  off  completely  just  before  the  film  run  to  eliminate  convection 
currents.  Lights  were  switched  on  manually  just  before  the  camera 
and  the  photographic  run  was  carried  out  at  the  preselected  speed. 

DISCUSSION  OF  RESULTS 

Observations  and  Comparison  with  Existing  Correlations 

Sixty-one  100-foot  rolls  of  high  speed  film,  taken  at  3,000  and 
5,000  frames  per  second,  were  analysed  frame  by  frame  using  a  Vanguard, 
Model  C-llD,  motion  picture  analyser.  Bubble  diameters  at  departure 
were  measured  at  the  instant  of  their  break-off  from  the  heating  surface. 
These  measurements  were  carried  out  for  either  all  bubbles  appearing  on 
a  film  or  a  limited  number  of  bubbles  considered  representative  of  the 
film  if  there  was  ?  large  number  of  bubbles  on  the  film.  For  each  high 
speed  film  taken,  the  pressure  at  the  heating  surface  (p),  the  heat  flux 
on  the  outer  surface  of  the  heating  element  (q),  the  average  outer  wall 
to  saturation  temperature  difference  (AT),  the  number  of  bubbles  measured, 
and  the  maximum,  minimum,  and  arithmetic  average  of  all  bubble  diameters 
measured  are  listed  separately  for  each  heating  element  in  Tables  1 
through  7.  The  maximum  and  minimum  bubble  departure  diameter  values 
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are  included  to  emphasize  variations  caused  by  factors  other  than 
those  considered  in  this  investigation.  Bubble  diameters  at  departure 
as  predicted  by  a  number  of  correlations  are  also  listed  in  these 
tables  to  permit  comparison  with  the  experimental  data.  In  tabulated 
order,  these  correlations  are  given  below.  Detailed  assumptions  upon 
which  derivation  of  these  correlations  are  based  were  included  earlier 
in  a  literature  survey  by  Garg.* 


For  cylindrical  heating  surfaces,  bubble  diameters  at  departure 
in  saturated  pool  boiling  as  predicted  by  Owens^,  Ivey  and  Morris^,  and 
Garg^,  respectively,  are  given  by 
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(3) 


where  D  =  bubble  diameter  at  departure  from  a  cylindrical 
heating  surface  (ft) 

dQ  =  outer  diameter  of  the  heating  element  (ft) 

C,  =  specific  heat  of  liquid  (Btu/lb  °F) 

■t  m 

3 

p  =  density  of  liquid  (lb  /ft  ) 

m 

3 

p  =  density  of  vapor  (lb  /ft  ) 
v  '  m  ' 

h^^  =  latent  heat  of  evaporation  (Btu/lb^) 

k  =  thermal  conductivity  of  liquid  (Btu/hr-f t-°F) 

C  =  experimental  constant  for  cylindrical  heating  surfaces 
which  relates  the  thickness  of  the  superheated  layer  in 
nucleate  boiling  where  a  bubble  envelops  the  heating 
surface  to  its  value  assuming  pure  conduction 
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For  cylindrical  heating  surfaces  where  the  departure  diameter  of 
bubbles  was  less  than  four  times  the  diameter  of  the  heating  surface, 
the  departure  diameter  of  bubbles  listed  in  the  tables  is  that  pre¬ 
dicted  by  the  following  equation  by  Garg4  for  a  flat  surface  instead 
of  Equation  (3),  as  shown  in  Tables  6  and  7. 


3  p^C^k ( AT ) 

4qP  h. 
v  fg 


(4) 


where  Df  =  maximum  bubble  diameter  at  departure  from  a  flat  surface 
(ft),  and  Cf  =  experimental  constant  which  relates  the  thickness  of 
the  superheated  layer  in  nucleate  boiling  to  its  value  assuming  pure 
conduction  for  flat  surfaces,  and  for  cylindrical  surfaces  where  the 
bubble  diameter  at  departure  is  less  than  four  times  the  diameter  of 
the  heating  surface. 

The  measured  average  bubble  diameters  at  departure  from  cylindrical 
heating  surfaces  where  the  bubble  diameters  exceeded  the  heating  element 
diameter  by  a  factor  of  four  were  converted  to  values  which  might  have 
been  obtained  from  flat  heating  surfaces  under  the  same  experimental 
conditions  using  an  equation  previously  developed  by  Garg4; 
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where  De  =  equivalent  maximum  bubble  diameter  at  departure  from 
a  flat  surface  (ft). 

The  value  of  Cc  was  determined  experimentally  and  found  to  be  2.0 
during  experiments  in  saturated  nucleate  pool  boiling  at  sub-atmospheric 
pressures  from  a  0 . 125- inch-OD  tube^.  The  values  of  both  Cc  and  C£  were 
assumed  to  be  2.0  for  the  purpose  of  comparisons  with  the  predicted 
values.  Finally,  the  last  column  in  the  tables  lists  bubble  diameters 
at  departure  from  flat  surfaces  as  predicted  by  Patten  and  Garg^,  using 
the  following  equation: 
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where 


o 

CT 


g 


critical  pressure  (lb^/in  abs) 

2 

conversion  factor  (lb  -ft/lb, -hr  ) 

m  f 

surface  tension  of  liquid  (lb^/ft) 

2 

acceleration  due  to  gravity  (ft/hr  ). 


The  maximum,  minimum  and  average  values  of  bubble  diameter  at 
departure  are  plotted  against  the  pressure  at  the  heating  surface, 
separately  for  each  heating  element  in  Figures  6  through  12  and  are 
compared  with  values  predicted  by  Equation  (3).  The  equivalent  bubble 
diameters  at  departure  from  a  flat  surface  as  obtained  by  modifying 
experimental  data  using  Equation  (5)  are  compared  with  the  predicted 
values  using  Equation  (6)  in  Tables  1  through  7.  Before  discussing  the 
difference  between  the  predicted  and  measured  values,  the  process  of 
bubble  formation,  growth  and  departure  may  be  described  which  might 
help  toward  explaining  the  scatter  in  the  data  and  their  observed 
deviation  from  the  predicted  values. 


a.  A  bubble  grows  by  absorbing  the  superheat  energy  from  the 
boundary  layer.  If  after  absorbing  the  available  superheat  energy  a 
bubble  has  not  been  able  to  reach  a  large  enough  size  to  overcome  the 
surface  tension  and  viscous  forces,  it  will  stay  attached  to  the 
heating  surface  and  continue  to  grow  by  evaporation  at  the  solid- 
liquid-vapor-  interface  ,  by  conglomeration  with  bubbles  growing  in  the 
vicinity,  or  both. 

b.  If  the  superheat  energy  in  the  boundary  layer  at  bubble 
initiation  is  large  enough  to  cause  this  bubble,  during  its  growth,  to 
exceed  the  size  necessary  to  overcome  the  surface  tension  and  viscous 
forces,  the  bubble  may  leave  the  heating  surface  even  though  its  growth 
has  not  yet  been  completed. 

c.  The  growth  rate  of  a  bubble  depends  upon  the  superheat  energy 
available  in  the  layer,  thermodynamic  and  transport  properties  of  the 
fluid,  conglomeration  with  other  bubbles  which  happen  to  be  in  the 
vicinity,  and  evaporation  at  the  interface  if  the  bubble  lias  not  attain¬ 
ed  a  large  enough  ,ize  to  leave  the  heating  surface. 

d.  During  growth,  a  bubble  begins  to  move  upward  as  soon  as  the 
buoyancy  forces  exceed  the  surface  tension  forces.  At  this  time,  the 
bubble  enters  into  a  dynamic  balance  with  the  buoyancy,  viscous  and 
surface  tension  forces.  Before  departure  from  the  heating  surface, 
the  bubble  center  has  to  travel  a  sufficient  distance  away  from  the 
heating  surface  to  clear  the  bottom  of  the  bubble  while  still  attached 
to  the  heating  surface. 
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e.  The  total  distance  travelled  by  a  bubble  in  a  given  time  from 
its  initiation  will  depend  upon  the  bubble  growth  rate  and  position 

of  the  bubble  on  the  heating  surface  at  the  time  of  initiation,  changes 
in  growth  rate  with  time  due  to  changes  in  the  residual  liquid  superheat 
energy  around  the  bubble,  bubble  activity  in  the  vicinity  and  conglomer¬ 
ation  of  other  bubbles  into  it. 

f.  The  superheat  energy  in  the  boundary  layer  depends  upon  the 
waiting  period  and,  to  some  extent,  upon  the  bubble  activity  in  its 
vicinity.  The  waiting  period  depends  strongly  upon  the  size  and  shape 
of  surface  cavities,  as  has  been  observed  by  a  number  of  investigators. 

Because  of  the  many  interdependent  factors,  especially  item  (f) 
above,  it  may  be  easily  seen  why  bubbles  vary  so  much  in  size  even  when 
externally  controllable  parameters  like  pressure,  heat  flux,  degree 
of  subcooling,  etc.,  are  held  constant. 

Observation  (a)  was  noted  with  smaller  diameter  wires.  Bubbles 
from  certain  nucleation  sites  were  found  to  stay  attached  to  the 
heating  surface  for  a  very  long  time,  and  continue  to  increase  in  size. 
For  example,  bubbles  growing  on  the  0.001-inch  diameter  wire  at  atmos¬ 
pheric  pressure  were  observed  to  stay  attached  to  the  wire  for  as  long 
as  233  milliseconds.  Such  long  attachment  times  were  experienced  by 
bubbles  only  if  they  did  not  conglomerate  with  other  bubbles.  At 
nucleation  sites  where  bubbles  conglomerated  during  their  attachment 
to  the  heating  surface,  the  attachment  time  was  much  smaller.  For 
the  larger  diameter  heating  surfaces,  the  bubble  diameters  were  much 
larger  and  the  attachment  time  much  smaller  than  the  smaller  diameter 
heating  surfaces.  For  the  largest  diameter  heating  surface,  it  appeared 
that  observation  (b)  may  well  apply. 

The  somewhat  large  deviation  of  bubble  diameters  from  values  pre¬ 
dicted  by  Equation  (3)  shown  in  Figures  6  and  7  may  be  explained  by 
the  fact  that  during  derivation  of  the  equation,  it  was  assumed  in  the 
interest  of  simplicity  that  the  thickness  of  the  superheated  boundary 
layer  in  boiling  is  small  compared  to  the  diameter  of  the  heating 
surface.  This  assumption  becomes  less  valid  as  the  lowest  diameter  of 
the  heating  surface  of  0.001  inch  is  approached.  Furthermore,  it  may 
be  noted  that  this  equation  was  obtained  from  a  simple  energy  balance 
and  as  such  does  not  consider  factors  to  satisfy  bubble  departure 
requirements.  Therefore,  for  small  diameter  wires  at  higher  pressures, 
the  superheat  energy  in  the  boundary  layer  may  not  be  enough  to  produce 
large  enough  bubbles,  resulting  in  their  staying  attached  to  the  heat¬ 
ing  surface  until  they  can  grow  to  a  sufficient  size  prior  to  departure, 
as  described  in  observation  (a).  For  tubes  of  0.0506  and  0.100- inch  OD, 
the  exact  opposite  may  well  have  begun  to  hold  true,  as  described  in 
observation  (b),  as  the  lowest  pressure  is  approached.  This  observation 
may  offer  an  explanation  for  measured  bubble  diameters  in  some  cases 
to  be  somewhat  smaller  than  the  values  predicted  by  Equations  (3)  and 
(4). 
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As  Che  system  pressure  is  increased  above  Che  atmospheric  pressure, 
it  is  easy  to  envision  that  bubble  growth  rates  will  be  reduced  sig¬ 
nificantly  because  of  two  factors;  lower  wall-to-saturation  temperature 
drops  in  boiling  which  will  limit  the  superheat  energy  in  the  boundary 
layer,  and  higher  vapor  densities  causing  a  smaller  final  size  for  the 
same  superheat  energy.  Under  such  conditions,  it  is  highly  probable 
that  the  growth  due  to  energy  in  the  boundary  layer  plays  a  progressively 
lesser  part  in  determining  the  final  size  at  departure.  The  well  known 
Fritz**  equation  may  well  begin  to  become  important  in  determining 
bubble  sizes  at  departure  at  the  higher  pressures.  However,  since  the 
Fritz  correlation  does  not  make  allowance  for  either  the  heat  transfer 
which  causes  bubbles  to  grow  large  at  the  low  pressures  or  the  dynamic 
nature  of  the  boiling  process  which  might  cause  bubbles  to  depart 
before  the  available  superheat  energy  in  the  boundary  layer  has  been 
completely  absorbed  by  the  bubble,  it  is  anticipated  that  bubble 
diameters  at  departure  may  well  be  above  or  below  the  values  predicted 
by  the  Fritz  correlation.  A  combination  of  the  Fritz  correlation  with 
consideration  of  the  superheat  energy  in  the  boundary  layer  may  well 
offer  a  better  correlation. 

Equations  (1)  to  (5)  are  based  only  upon  an  energy  balance  between 
the  superheated  liquid  layer  and  the  bubble,  and  they  do  not  include 

the  requirement  that  a  bubble  has  to  continue  to  grow  to  a  large 

enough  size  to  be  able  to  overcome  the  surface  tension  and  viscous 
forces  prior  to  its  departure  from  the  heating  surface.  An  examination 

of  Figures  6  to  12  and  Tables  1  to  7  shows  that,  despite  this  limitation 

and  the  unavoidable  scatter  of  data,  these  equations  are  much  more 
successful  in  predicting  bubble  diameters  than  the  Fritz  correlation. 

This  suggests  that  the  superheat  energy  in  the  boundary  layer  in  these 
experiments  must  have  been  a  much  more  significant  factor  in  determining 
bubble  diameters  at  departure  compared  to  the  effects  of  buoyancy  and 
surface  tension  forces,  at  least  at  the  lower  pressures.  Some  adjustment 
in  these  equations  to  include  the  presence  and  the  requirement  of  a  suffi¬ 
cient  buoyancy  force  is  needed  for  predictions  at  the  higher  pressures. 

The  two  factors  discussed  above  viz,  (1)  the  assumption  of  thick¬ 
ness  of  the  boundary  layer  being  small  compared  to  the  diameter  of  the 
heating  surface  and,  (2)  the  bubble  departure  requirements,  may  explain 
the  observed  increase  in  deviation  between  the  predicted  and  the 
measured  bubble  diameters  as  the  heating  element  diameter  decreases. 

This  was  especially  true  at  the  higher  pressures  where  bubbles  stayed 
attached  to  the  heating  surface  and  kept  growing  until  they  reached  a 
large  enough  size  to  permit  departure  from  the  heating  surface. 

The  systematic  effects  of  system  parameters  considered  as 
variables  in  this  experimental  investigation  will  now  be  discussed. 


Effect  of  Diameter  of  the  Heating  Surface 

Bubble  diameters  at  departure  from  a  heating  surface  at  a  partic¬ 
ular  pressure  were  averaged  for  all  heat  fluxes  and  plotted  against 
diameter  of  the  heating  element,  separately  for  each  pressure  in 
Figure  13.  This  was  done  because,  as  will  be  shown  later,  the  effect 
of  heat  flux  upon  bubble  diameter  was  undetectably  small  compared  to 
variations  caused  by  factors  other  than  those  considered  in  this 
investigation.  Also  shown  in  the  figure  for  comparison  is  the  dependence 
of  bubble  diameter  upon  heating  element  diameter  as  predicted  by 
Equation  (3),  viz., 

D  «dJ  . 
c  o 

If  allowance  is  made  for  the  observation  discussed  earlier  that  bubbles 
from  smaller  diameter  wires  were  much  larger  because  of  the  departure 
requirements,  the  data  approximately  support  the  dependence  predicted 
by  Equation  (3). 

Effect  of  Pressure 

Average  bubble  diameters  at  departure  are  plotted  against  pressure 
separately  for  each  heating  element  in  Figure  14.  As  may  be  seen,  the 
bubble  diameter  decreases  w:.th  an  increase  in  the  system  pressure. 
Although  there  is  no  pressure  term  in  Equation  (3),  its  presence  is 
well  accounted  for  through  its  effect  upon  the  properties  of  liquid 
and  vapor.  A  more  direct  inclusion  of  the  pressure  term  in  a  correlation 
was  carried  out  in  Equation  (6). 

As  shown  later,  a  correlation  of  all  data  was  carried  out  by 
combining  Equation  (3)  with  the  Fritz  correlation  to  include  the 
effect  of  buoyancy  and  surface  tension  forces. 

Effect  of  Heat  Flux 

To  determine  the  systematic  effect  of  heat  flux  upon  the  bubble 
diameter  at  departure,  the  ratio  of  average  bubble  diameters  at 
departure  for  the  highest  and  lowest  heat  flux  values  was  compared 
with  the  ratio  of  the  highest  and  the  lowest  heat  fluxes.  These 
ratios  are  tabulated  for  all  wires  and  tubes  at  all  system  pressures 
in  Table  8.  Also  given  in  this  table  for  comparison  is  the  highest 
value  of  the  ratio  of  maximum  and  minimum  bubble  diameters  at  a  given 
heat  flux.  In  Table  8,  therefore 

Average  (D  )  for  the  highest  q 

B  s  9 
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It  may  be  seen  from  the  values  of  ratio  B3  in  Table  8  that,  for 
a  given  set  of  experimental  conditions  (pressure  at  the  heating  sur¬ 
face,  diameter  of  the  heating  element,  surface  orientation,  average 
heat  flux,  height  of  water  above  the  heating  surface,  and  the  degree 
of  subcooling  which  was  zero  in  the  work  reported  here),  the  variation 
in  bubble  diameter  at  departure  was  as  high  as  X4.10  for  the  0.100-inch- 
OD  tube  at  3-0  psia.  It  may  also  be  noted  that  with  this  one  exception 
in  the  value  of  ratio  B3  of  4.10  for  the  0 . 100- inch-OD  tube  at  3.0 
psia,  the  ratio  B3  decreases  as  the  system  pressure  is  increased. 

This  anomaly  was  apparently  caused  by  large  waiting  periods  on  the 
entire  tube  surface  followed  by  explosive  growth  of  large  bubbles 
at  a  specific  location  on  the  heating  surface.  At  departure,  these 
large  bubbles  left  vapor  tails  on  the  tube  surface  which  grew  into 
trails  of  smaller  bubbles  from  these  locations.  The  former,  parent 
bubbles  were  perhaps  the  result  of  nucleation  of  a  surface  cavity 
which  required  a  higher  wall  superheat,  thereby  a  larger  superheat 
energy  at  bubble  initiation,  whereas  the  streams  of  bubbles  following 
these  large  bubbles  were  probably  caused  by  evaporation  at  the 
boundary  of  the  pre-existing  gaseous  phase.  These  latter  bubbles  had 
very  little  waiting  period  and  smaller  attachment  time  to  the  heating 
surface  compared  to  the  parent  bubbles.  It  was  perhaps  this  difference 
in  nucleation  of  a  surface  cavity  in  one  case  and  evaporation  into 
the  pre-existing  gaseous  phase  in  the  other  which  caused  sucli  large 
differences  in  their  departure  sizes. 

Variations  in  bubble  size  at  departure  could  be  caused  by  a 
combination  of  several  factors.  Besides  the  waiting  period  and  the 
attachment  time  of  the  bubble  to  the  heating  surface  discussed  above, 
other  important  parameters  which  may  affect  maximum  bubble  size  are: 
variations  in  the  size,  shape  and  density  of  cavities  in  the  heating 
surface,  the  bubble  growth  rate  and  the  dynamics  of  bubbles  growing 
in  the  vicinity. 

In  order  to  evaluate  the  effect  of  heat  flux  upon  bubble  diameter 
at  departure,  two  factors  will  now  be  examined.  First,  the  ratio  B] 
which  was  found  to  be  below  1.00  for  11  of  the  23  test  conditions 
indicating  that  an  increase  in  heat  flux  reduces  the  bubble  diameter 
at  departure,  and  above  1.00  for  the  other  12  test  conditions  indicating 
exactly  the  opposite.  Second,  with  one  exception  the  variation  of 
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ratio  Bi  above  and  beLow  1.00  was  found  to  be  much  smaller  than  variation 
in  bubble  diameter  at  departure  indicated  by  the  ratio  Bj.  An  examina¬ 
tion  of  these  two  factors  indicates  little  if  any  effect  of  heat  flux 
upon  the  maximum  bubble  diameter  at  departure.  It  may  well  be  that 
the  heat  flux  does  have  some  systematic  effect  upon  the  bubble  diameter; 
however,  in  view  of  the  large  values  of  ratio  B3,  its  determination 
will  have  to  wait  until  the  effects  of  parameters  other  than  those 
considered  in  this  work  upon  the  bubble  diameter  at  departure  are 
systematically  determined.  For  this  analysis,  therefore,  the  effect 
of  heat  flux  will  be  considered  to  be  negligible.  This  observation  is 
in  agreement  with  the  work  of  Gaertner  and  Westwater?  and  Hatton  and 
Hall.8 

Correlation  of  Experimental  Data 

If  the  above  observation  of  negligible  effect  of  heat  flux  on 
bubble  diameter  at  departure  is  assumed  to  be  reasonably  accurate 
over  the  entire  discrete  bubble  region  at  a  given  pressure,  and  since 
only  one  of  heat  flux  or  wall  to  bulk  temperature  difference  is  inde¬ 
pendently  variable  at  a  constant  pressure,  the  heat  flux  may  be  assumed 
to  be  proportional  to  the  square  of  the  wall  to  bulk  temperature 
difference,  or 


q  =  Cx  (4 T)2  .  (10) 


If  Equation  (10)  can  be  combined  with  Equations  (3)  and  (4),  the  need 
for  experimental  measurements  in  predicting  bubble  diameters  at 
departure  may  be  eliminated.  However,  some  allowance  has  to  be  made 
for  this  constant  of  proportionality  which,  although  assumed  a  constant 
over  the  entire  discrete  bubble  region  at  a  given  pressure,  is  not  a 
constant  for  all  pressures  This  adjustment  may  be  carried  out  as 
follows.  Since  all  the  liquid  and  vapor  properties  used  in  Equations 
(3)  and  (4)  vary  with  pressure,  the  variation  in  the  value  of  C3  may 
be  accounted  for  by  changing  the  index  of  all  the  properties  combined. 
Therefore,  Equations  (3)  and  (4)  modified  by  Equation  (10)  may  be 
written  as 


D2  =  3d  C 
c  o  c 


(11) 
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and 


D 


f 


K 


(12) 


where  and  C2  are  experimental  constants.  Combining  Equations  (11) 
and  (12)  leads  to 


D 


f 


(13) 


which  is  of  the  same  form  as  Equation  (5)  and  may  be  used  to  convert 
bubble  diameters  at  departure  from  cylindrical  surfaces  into  equiva¬ 
lent  bubble  diameters  at  departure  from  flat  surfaces. 

To  include  the  effects  of  buoyancy  and  surface  tension  forces 
which  determine  the  bubble  departure  requirements,  and  assuming  the 
contact  angle  to  be  constant  at  all  pressures^,  Equation  (12)  may  be 
modified  by  the  Fritz  expression  to  give 


D 


f 


2  c  c 

4  3  Cf 


~  w  1 2  r  s°g 

LcipvhfgJ  Lg(vpv>J 


(14) 


where 


is  a  constant. 


Equation  (14)  may  be  simplied  to 


r  w 

r  s°a  1 

KhtgJ 

_g<pt-pv>_ 

(13) 


where  C,  is  a  dimensional  constant. 

4 

All  experimentally  obtained  bubble  diameters  at  departure  for 
water  from  cylindrical  wires  and  tubes  where  the  bubble  diameter 
exceeded  the  diameter  of  the  heating  surface  by  a  factor  of  four  were 
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modified  using  Equation  (5)  and  compared  with  Equation  (15)  in  Figure 
15.  Experimental  data  from  other  authors^*®"^  for  water  taken  in 
the  discrete  bubble  region  were  also  inserted  in  the  figure,  and  a 
suitable  value  for  C2  was  found  to  be  1.05.  The  value  of  constant  C4 
will  depend  upon  the  units  used  in  Equation  (15).  The  large  spread  of 
data  in  Figure  15  is  not  unexpected  because, 

a.  the  heat  flux  is  not  always  proportional  Lo  the  square  of 
wall  superheat, 

b.  the  size,  shape  and  distribution  of  cavities  on  the  heating 
surface  have  a  large  effect  upon  the  superheat  energy  in  the  boundary 
layer  at  the  time  of  bubble  initiation  through  their  effects  upon  the 
waiting  period, 

c.  the  thickness  of  the  boundary  layer  and  the  temperature  grad¬ 
ient  in  it  in  boiling  cannot  be  accurately  predicted  by  the  use  of  a 
simple  conduction  equation, 

d.  heat  flux  and  wall  to  saturation  temperature  difference  in 
boiling  are  not  static  quantities,  and  they  vary  greatly  in  time  and 
location  on  the  heating  surface,  and 

e.  bubble  activity  in  the  vicinity  of  a  nucleation  site  has  a 
significant  effect  upon  its  nucleation  characteristics  and  upon 
bubble  growth  and  departure  characteristics. 

Since  not  one  of  the  above  five  can  be  quantified  in  a  meaningful 
manner,  accurate  prediction  of  bubble  diameter  is  extremely  difficult. 
Nevertheless,  attempts  to  at  least  determine  the  gross  effects  of  the 
thermodynamic  and  transport  properties  of  the  liquid  and  vapor  have  to 
be  made.  Hopefully,  Equation  (15)  is  a  step  in  the  right  direction. 

The  next  step  is  to  experimentally  investigate  the  gross  effects 
of  roughness  of  the  heating  surface  and  small  variations  in  bulk  temp¬ 
eratures  from  the  saturation  temperature  upon  the  bubble  diameters  at 
departure . 

CONCLUSIONS 

A  high  speed  photographic  study  of  saturated  nucleate  pool  boiling 
of  water  was  carried  out  from  horizontally  mounted  wires  and  tubes  of 
diameter  between  0.001  and  0.100  inch  at  subatmospheric  pressures. 

Based  upon  the  frame  by  frame  analysis  of  films,  the  following  conclu¬ 
sions  were  reached. 

a.  The  bubble  diameter  at  departure  may  be  determined  from  an 
energy  balance  in  the  superheated  boundary  layer,  with  some  considera¬ 
tion  for  the  buoyancy  and  surface  tension  forces  which  control  bubble 
departure . 
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b.  The  diameter  of  a  bubble  at  departure  is  approximately  pro¬ 
portional  to  the  square  root  of  the  heating  surface  diameter  for  tests 
where  the  bubble  diameter  at  departure  exceeds  the  heating  surface 
diameter  by  a  factor  of  four. 

c.  An  increase  in  the  system  pressure  decreases  the  diameter  of 
the  bubble  at  departure. 

d.  There  is  no  systematic  effect  of  heat  flux  upon  the  diameter 
of  bubbles  at  departure  in  the  discrete  bubble  region. 

e.  An  empirical  equation  is  proposed  which  successfully  corre¬ 
lates  all  the  available  data  on  bubble  diameters  at  departure  in 
saturated  nucleate  pool  boiling  of  water  at  all  pressures. 

RECOMMENDATIONS 

A  systematic  experimental  investigation  should  be  carried  out 
to  determine  the  effects  of  surface  roughness  and  small  variations  in 
bulk  temperature  from  the  saturation  temperature  upon  the  bubble  dia¬ 
meters  at  departure  in  pool  boiling. 


Table  3.  Data  for  Saturated  Pool  Boiling  of  Water 
for  a  Wire  of  Diameter  0.005  Inch 
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Table  5.  Data  for  Saturated  Pool  Boiling  of  Water 
for  a  Wire  of  Diameter  0.026  Inch 


Table  6.  Data  for  Saturated  Pool  Boiling  of  Water 
for  a  Tube  of  Diameter  0.0506  Inch 
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Obtained  from  Equation  (4)  because  the  bubble  diameter  was  less  than  4  times  the  diameter 
of  the  heating  surface. 

Same  as  experimentally  measured  diameter  since  it  was  less  than  4  times  the  diameter  of  the 
heating  surface. 


Table  7.  Data  for  Saturated  Pool  Boiling  of  Water 
for  a  Tube  of  Diameter  0.100  Inch 
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Obtained  from  Equation  (4)  because  the  bubble  diameter  was  less  than  4  times  the  diameter  of 
the  heating  surface. 

Same  as  experimentally  measured  diameter  since  it  was  less  than  4  times  the  diameter  of  the 
heating  surface. 
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Table  8. 


Data  for  the  Determination  of  the  Effect  of 
Heat  Flux  Upon  Bubble  Diameter  at  Departure 
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*  Data  not  available. 
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1  kw  diffused 


600  w  incandescent 
lights 


Fig.  2.  Arrangement  of  lights  and  camera  for  high 
speed  photography. 


220V 

AC 


C  Condenser 
D  Drain  valve 
DCP  D.C.  power  supply 
H  Auxiliary  heater 
LV  Controlled  leak 
needle  valve 
M  Manometer 
P  Potentiometer 
R  Resistance 
S  Selector  switch 
T  Thermocouple 
VD  Voltage  divider 
VP  Vacuum  pump 


220VAC 


Fig.  3.  Schematic  of  the  experimental  arrangement. 
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Fig.  4.  Photograph  of  the  experimental  arrangement. 
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Fig.  5.  Photograph  of  the  experimental  arrangement 
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predicted  by  Equation  (3). 
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p,  psia 


Fig.  7.  Maximum,  minimum  and  average  values  of  bubble 

diameters  at  departure  from  0.002-inch-diameter 
wire  at  various  pressures  and  heat  fluxes. 
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Fig.  10.  Maximum,  minitnuu  and  average  values  of  bubble 

diameters  at  departure  from  0.026-inch-diameter 
wire  at  various  pressures  and  heat  fluxes. 
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p,  psia 

Fig.  14.  Effect  of  pressure  upon  the  average  bubble 
diameter  at  departure  in  saturated  nucleate 
pool  boiling  of  water. 
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Symbol  Ref. 
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10 

□ 

11 

o 

12 

o- 

13 

(Equilibrium  bubbles) 

13 

(Coalesence  bubbles) 

#- 

13 

(Boiling  bubbles) 

X 

14 

+ 

15 

■ 

This  work 

10 


10 


10 

P.  psia 


10' 


Fig.  15.  Comparison  of  modified  maximum  bubble  diameters  with  the 
correlation  for  nucleate  pool  boiling  of  saturated  water. 
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APPENDIX  A 


MEASUREMENT  OF  DIAMETER  OF  THE  VARIOUS  HEATING 
ELEMENTS  USED  IN  THE  TEST  PROGRAM 

0.001-,  0.002-  and  0.005-inch-nominal-diameter  wires 

The  diameters  of  these  wires  were  considered  too  small  to  be 
measured  accurately  either  by  a  micrometer  or  by  weighing  known 
lengths  using  an  accurate  balance.  The  diameters  of  these  wires  were 
therefore,  determined  using  the  resistivity  method.  The  values  of 
resistivity  and  temperature  coefficient  of  resistivity  were  obtained 
from  measurements  carried  out  on  the  0.02595-inch-diameter  wire.  The 
details  of  resistivity  measurements  are  given  in  Appendix  B.  From 
these  measurements,  the  diameters  of  the  0.001-,  0.002-  and  0.005- 
inch-nominal-diameter  wires  were  found  to  be  0.00107,  0.00201  and 
0.00500  inch,  respectively. 

0.010- and  0.025-inch-noMinal-diameter  wires 

The  diameters  of  these  wires  were  considered  too  small  to  be 
measured  accurately  by  a  micrometer.  Since  these  wires  were  heavy 
enough  to  permit  determination  of  their  diameters  by  accurately  weigh 
ing  known  lengths,  this  method  was  preferred  over  the  less  accurate 
resistivity  method.  The  lengths  and  weights  of  approximately  nine 
inch  long  samples  were  measured  accurately.  Assuming  the  density  of 
platinum  to  be  1335  pounds  per  cubic  foot,  these  diameters  were 
calculated  to  be  0.01005  and  0.02595  inch,  respectively. 

0.05-  and  0 . 10-inch-nominal-0D  tubes 

Accurate  measurement  of  the  outer  diameters  of  Lhe  0.05-  and 
0. 10-inch-nominal-OD  tubes  was  carried  out  using  a  micrometer  capable 
of  reading  0.0001  inch.  Average  OD  of  the  two  tubes  was  found  to  be 
0.0506  and  0.1000  inch,  respectively.  The  average  inner  diameter  of 
these  tubes  was  determined  indirectly  by  accurately  weighing  known 
lengths,  and  found  to  be  0.03725  and  0.09116  inch,  respectively. 
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APPENDIX  15 

MEASUREMENT  OF  RESISTIVITY  AND  TEMPERATURE  COEFFICIENT 
OF  RESISTIVITY  OF  THE  VARIOUS  WIRES  AND  TUBES 

Approximately  six  inch  lengths  of  various  wires  and  tubes  were 
mounted  horizontally  between  current  carrying  conductors  and  the 
lengths  measured  accurately  by  a  vernier  caliber  capable  of  measuring 
it  to  an  accuracy  of  0.001  inch.  The  assembly  was  lowered  into  a 
tank  of  demineralized  water  at  room  temperature  and  the  wire  or  tube 
placed  as  resistance  R  in  the  circuit  shown  below. 


where  R  -  standard  resistance  (ohm) 
s 

R  =  200  Q  (approx.)  for  wires  of  nominal  diameter 
between  0.001  and  0.010  incli 

-  27  fl  (approx.)  for  wire  of  0.025-inch  nominal 
diameter  and  tubes  of  0.05-  and  0.10-inch 
nominal  OD 

R  =  resistance  of  the  heating  element  under 
investigation  at  temperature  T  (ohm) 

e  ,  e  =  voltage  drop  across  R  and  R  ,  respectively  (volt). 

S  X  S  X 

The  resistance  R  was  used  simply  to  reduce  the  power  input  to 
resistances  Rx  and  Rs.  The  temperature  of  water  in  the  tank  was 
maintained  at  either  room  temperature  or  at  the  boiling  point  at 
atmospheric  pressure,  and  measured  accurately  in  each  case  at  the 
wire  or  tube  level  by  a  calibrated  thermocouple . 

During  calibration  of  0.001- ,  0.002-,  0.005-  and  0.010-inch- 
nominal -diameter  wires,  the  emf  from  the  dc  source  was  adjusted  to 
limit  the  current  flow  to  between  20  and  25  milliamperes .  Similarly, 
for  the  0.025-inch-nominal-diameter  wire  and  the  0.005-  and  0.10- 
inch-nominal-OD  tubes,  the  current  flow  was  limited  to  between  140 


43 


and  190  milliamperes .  At  these  low  current  levels, 
rise  in  the  wires  or  tubes  was  limited  to  less  than 
steady  state  conditions  were  reached,  the  values  of 
temperature  were  recorded. 


the  temperature 

0 . 1°F .  After 

e  ,  e  and  water 
s  x 


For  wires  of  nominal  diameter  0.001,  0.002  and  0.005  inch,  the 
accurate  values  of  the  diameter  were  obtained  from: 


o 


4032^es 

TTR  e 
s  x 


1  + 


a(Tf-32) 


(B-l) 


where  d 

o 


l 


a 


T 


f 


outer  diameter  of  the  heating  element  (ft) 

resistivity  of  the  heating  element  at  32°F  (ohm-ft) 

length  of  the  heating  element  (ft) 
temperature  coefficient  of  resistivity  (F  *) 
temperature  of  water  (F) 


The  values  of  0^2  an^  a  f°r  Equation  (B-l)  were  obtained  from 
tests  carried  out  with  the  0.02595-inch-diameter  wire.  These  values 
were  $32  -  3.38  x  10"^  (ohm-ft)  and  a  =  0.0020  (f~^). 

The  resistivity  and  temperature  coefficient  of  resistivity  values 
for  wires  of  0.01005  and  0.02595  inch  diameter  were  obtained  by  solving 
the  simultaneous  Equations  (B-2)  and  (B-3),  and  those  for  the  0.0506- 
and  0. 100-inch-OD  tubes  were  obtained  by  solving  Equations  (B-4)  and 
(B-5),  shown  below,  all  using  measured  values  of  eg  and  ex  at  two 
temperatures.  For  wires, 


!¥*!  t -  ^{i+«<w32>}  • 


(B-2) 


R  e 
s  x 


%2l 


B  ’  ht  TTd 


1  +  a<Tfht~32)}  : 


(B-3) 
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r 


for  tubes. 


(B-4) 


(B-5) 

'  o  1" 

where  d^  is  the  inner  diameter  of  tube  (ft),  subscript  rt  is  the  value 
at  room  temperature,  and  subscript  ht  is  the  value  when  boiling  at 
atmospheric  pressure. 

The  following  computed  values  were  used  in  all  subsequent 
calculations . 


O  1 


R  e 


s  ’  ht  TT(d  -d  .  ) 


40  2*”  i 

1-Xl1  +  a<Tfht-32)^ 


wire  or 
tube  OD 
(in) 

^32 

(ohm-f t) 

1 

&  Ui 

o 

0.00107 

3.38  x  10~7 

0.0020 

0.00201 

3.38  " 

0.0020 

0.00500 

3.38  " 

0.0020 

0.01005 

3.325  " 

0.0020 

0.02595 

3.38 

0.0020 

0.0506 

3.681  " 

0.002016 

0.1000 

3.715  " 

0.001874 
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APPENDIX  C 


DETERMINATION  OF  SURFACE  TEMPERATURE  OF 
THE  HEATING  ELEMENTS 

In  all  tests  reported  here,  the  outer  surface  temperature  of  heating 
elements,  Tw,  was  calculated  from  the  measured  values  of  the  current 
flowing  through  and  the  voltage  drop  across  them.  Spot  welding  a 
thermocouple  on  the  wire  or  tube  surface  was  considered  unsuitable  for 
this  investigation  because,  (a)  spot  welding  a  thermocouple  would  have 
significantly  changed  the  characteristics  of  the  location  in  boiling 
to  make  measurements  unrepresentative  of  the  entire  surface,  (b)  the 
wire  diameter  and  the  tube  wall  thicknesses  were  small,  (c)  only  a 
small  surface  area  of  the  thermoucouple  would  have  actually  come  in 
contact  with  the  heating  surface,  the  rest  being  in  contact  with  the 
liquid,  and  (d)  local  wall  temperature  fluctuations  caused  by  boiling 
would  have  made  measurement  of  the  average  heating  surface  temperature 
very  difficult.  Instead,  it  was  decided  to  determine  the  average  wall 
temperature  using  the  heating  surface  as  a  resistance  heater,  as 
described  below. 

For  a  given  set  of  experimental  conditions,  the  measured  values 
of  the  current  flowing  through  the  heating  element  and  the  voltage 
drop  across  it  gave  the  value  of  the  total  resistance  of  the  heating 
element  at  its  average  or  mean  temperature.  From  this  average  temp¬ 
erature  value,  the  outer  surface  temperature  was  calculated  as  follows. 

Wires 


Assuming  the  temperature  to  vary  with  radius  only,  the  heat  gen¬ 
erated  in  a  cylinder  of  radius  r  is  conducted  outwards,  and  may  be 
written : 


q'"  n r2l  =  -2nr£K 

n  dr 


(C-l) 


where  q"'  is  the  heat  generation  rate  in  the  heating  element  (Btu/hr-cu  ft), 
K  is  the  thermal  conductivity  of  the  heating  element  (Btu/hr-f t-F) ,  and 
T  is  the  temperature  (F).  If  T  =  Tc  at  r  =  0,  and  T  =  Tw  at  r  =  rQ, 
then  Equation  (C-l)  may  be  integrated  to  give: 


(C-2) 


46 


For  small  diameter  wires,  the  temperature  drop  from  the  center  of 
the  wire  to  outside  is  not  large.  Therefore,  in  this  derivation,  it 
was  assumed  that  the  effect  of  temperature  coefficients  of  electrical 
resistivity  and  thermal  conductivity  upon  the  heat  generation  rate  in 
the  heating  element  was  negligible.  The  heat  generation  rate  in  the 
wire  may  now  be  equated  to  the  heat  flux  on  its  surface  to  give: 


q'"nr  K>  =  q2TTr  t 
o  o 


or 


q 


/// 


From  Equations  (C— 2 )  and(C-4), 


T 

c 


9 


(C-3) 


(C-4) 


(C  -5) 


where  q 


3.413  V  1 

2TTr  £ 
o 


(C-6) 


V  =  voltage  drop  across  the  heating  element  (volt) 

I  =  current  flowing  through  the  heating  element  (amp). 


To  confirm  the  assumption  of  small  temperature  drop  within  the 
wire,  temperature  drop  across  the  largest  diameter  wire  used  in  these 
tests  was  calculated  at  the  maximum  experimental  heat  flux  of  54,600 
Btu/hr-sq  ft,  and  found  to  be  0.7°F.  Therefore,  the  assumption  of 
constant  heat  generation  rate  was  reasonable.  The  mean  temperature  of 
the  wire  may  now  be  related  to  its  outer  wall  temperature  as  follows. 
Assuming  a  linear  temperature  drop  across  the  wire  for  small  diameter 
wires,  the  wire  temperature  at  any  radius  r  may  be  written  as: 


r  -r 

T  =  T  +  — —  (T  -  T  )  , 

w  r  c  w'  ’ 
o 


and  the  mean  temperature  of  the  wire  may  be  written  as: 


(C-7) 


*7 


T  =■ 
m 


/o  r  -r  « 

2nrdr  <T  +  — —  (T  -  T  )> 
)  w  r  c  w'f 
o  1  o  ' 


TTr 


-  T  4  —  (T  -  T)  . 
w  3  c  w 


(C-8) 


The  mean  temperature,  Tm,  is  calculated  from  measured  values  of  current, 
voltage  drop,  and  length  and  diameter  of  the  wire  as: 


TTr  V  -  I <t„0l 

T  i32+  -° . 32 

Tm  ai0~C 


(C-9) 


From  Equations  (C-5),  (C-6),  (C-8)  and  (C-9),  one  obtains: 


T  = 
w 


nr2  V 
o 


i032i 


oil<t>32l 


4  32  - 


3.413  V  I 
12TTK*. 


(C-10) 


The  corrections  in  the  value  of  Tw  caused  by  the  temperature  co¬ 
efficients  of  thermal  conductivity  and  electrical  resistivity  were 
found  to  be  negligible. 

Tubes 


Assuming  linear  variations  in  thermal  conductivity  and  electrical 
resistivity  with  temperature,  the  temperature  drop  across  the  wall  of 
a  tube  may  be  written^  as: 


T.-T  = 


1  “  kiK/'i>2-i}L 


qr  (R  /R.) 
o  m 


6Ki 


i^  / 1  o ~ 1  i  1  (  3y44o'YTi40  > 

T2“ TT  l  I  i(l4aT.)  (14YT.)/+  ••• 


,  (C-ll) 


where  T^  =  inner  tube  wall  temperature  (F) 
r.  =  inner  tube  radius  (ft) 

l 
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R  =  mean  tube  resistance  (ohm) 
m 

-  tube  resistance  at  the  tube  inner  wall  temperature  (ohm) 

K.  =  thermal  conductivity  of  the  tube  at  inner  tube  wall 
temperature  (Btu/hr-f t-F) 

y  =  temperature  coefficient  of  thermal  conductivity  (F  ^). 
Simplifying  the  geometric  term,  Equation  (C-ll)  may  be  written  as: 


T.-T 
1  w 


"  ro  <yy 

Ki  {‘W2-1} 


‘V'.> 


q  r  (R  /R.)  .r  -r 

1  /-\  m  l  I  a 


3y4QryT.-fa 


+ 6ki  <iwii>  ‘i+vii>^ 


(C-12) 


To  determine  the  influences  of  the  temperature  coefficients  of 
electrical  resistivity  and  thermal  conductivity,  the  temperature  drop 
across  the  tube  wall  may  now  be  calculated  for  both  tubes.  For  the 
0. 100-inch-OD  tube  at  the  highest  value  of  heat  flux  of  55,250 
Btu/hr-sq  ft  at  1.0  psia  obtained  in  the  test  program,  the  tempprature 
drop  across  the  tube  wall  assuming  a  -  0  and  Y  =  0  was  found  to  be 
0.25°F.  Using  values  of  a  =  0.0020°F"^  and  y  =  0.0003°F-^,  the 
correction  term  was  found  to  be  negligibly  small.  Similarly,  the 
temperature  drop  across  the  wall  of  the  0 . 0506- inch-OD  tube  for  the 
highest  heat  flux  of  44,850  Btu/hr-sq  ft  at  1.0  psia  was  found  to  be 
0.32°F,  the  correction  term  involving  a  and  y  being  negligible. 


Since  only  the  mean  temperature  of  the  tube  wall  was  determined 
from  resistance  measurements,  values  of  (T^  -  Tw)  obtained  from 
Equation  (C-12)  had  to  be  modified.  The  mean  temperature  was  related 
to  the  tube  outer  wall  temperature  as  follows.  Referring  to  the 
tube  cross-section  shown  below,  the  temperature  at  any  radius  r  is 
given  by: 
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(C- 13) 


/  r  -r  \ 

T  =  T  4-  (T.  -  T  )  <  — -  >  , 

w  1  w  ( r  -r .  J 


o  1 


and  the  mean  temperature  of  the  tube,  Tm,  may  be  written  as: 


T  = 
m 


J  2nrdr  /T  +  (T.-T  )  /— - U 

ri  rw  1  w  \vri/> 


it  +  (t. 

w 

.  ~2  2, 

TT(r  -  r.) 
o  1 


T  +  (T. 
w  l 


T  ) 
w 


2  2 
r  4-  r .  r  -2  r . 
o _ 1  o _ 1 

2  2 

3  <ro  -  ri> 


(C- 14) 


The  values  of  the  terms  in  the  bracket  for  the  0.0506-inch-OD 
and  0. 100-inch-OD  tubes  were  found  to  be  0.47  and  0.49,  respectively. 
Therefore,  Equation  (C-14)  may  be  written  as: 

T  =  T  4  0.47  (T.  -  T  )  for  the  0 . 0506- inch-OD  tube  (C- 15) 
m  w  l  w 


and 


T  =  T  4  0.49  (T.-T  )  for  the  0. 1000-inch-OD  tube  .  (C-16) 

m  w  i  w 

Based  upon  the  measured  values  of  current  flowing  through  and  the 
voltage  drop  across  a  tube  of  length  l>,  the  mean  temperature  may  be 
written  as : 
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(C-17) 


n(r2-r2)  V  -  U  l 
-,0,  01  32 

T  -  32  +  - — — T - 

m  al®22<' 

Ignoring  Che  correction  term  due  to  temperature  coefficients  of 
electrical  resistivity  and  thermal  conductivity  and  assuming  R  =  R. , 
Equation  (C-12)  gives: 

T.-T  =  2.92  x  10  4  q/K.  for  the  0.0506-inch-OD  tube  (C-18) 


and 

-4 

T.-T  =  1.86  x  10  q/K.  for  the  0.  lOOO-inch-OD  tube  .  (C-19) 
l  w  1 

The  outer  wall  temperature  of  the  0.0506-inch-OD  tube  may  now  be 
obtained  by  combining  Equations  (C-15),  (C-17)  and  (C-18)  to  give: 


T  =  32  + 
w 


TT(r2-r2)  V  -  1 
o  i _ 32 

°^32^ 


1.37  X  10_4q 
k! 


(C-20) 


and  for  the  0 . 1000- inch-OD  tube,  Equations  (C-16),  (C-17)  and  (C-19) 
may  be  combined  to  give: 


T 

w 


32  + 


"<Vri>  V  •  U12l 


or 


1032^ 


9.11  x  10_5q 

T. 


(C-21) 


For  all  data  reported  here,  the 
determined  using  Equation  (C-10)  and 
0. lOOO-inch-OD  tubes  using  Equations 


surface  temperature  of  wires  was 
those  of  the  0.0506-inch-  and 
(C-20)  and  (C-21),  respectively. 
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NOMENCLATURE 


Surface  area  of  the  heating  element 
Ratio  defined  by  Equation  (7) 


ft 


Ratio  defined  by  Equation  (8) 


Ratio  defined  by  Equation  (9) 


Cj.Cj.Cj 


Cons  tants 


Dimensional  Constant 


Constant  for  cylindrical  heating  surfaces  which 
relates  the  thickness  of  the  superheated  layer 
in  nucleate  boiling  where  a  bubble  envelops  the 
heating  surface  to  its  value  assuming  pure 
conduction,  as  defined  by  Equation  (3) 

Constant  which  relates  the  thickness  of  the  super¬ 
heated  layer  in  nucleate  boiling  to  its  value 
assuming  pure  conduction  for  flat  surfaces  and 
for  cylindrical  surfaces  where  the  bubble 
diameter  at  departure  is  less  than  four  times 
the  diameter  of  the  heating  surface,  as  de¬ 
fined  by  Equation  (4) 


d. 

l 


Specific  heat  of  liquids 


Inner  diameter  of  tube 


Btu/lb  °F 
m 


ft 


Outer  diameter  of  the  heating  element 

Bubble  diameter  at  departure  from  a 
cylindrical  heating  surface 

Equivalent  maximum  bubble  diameter  at 
departure  from  a  flat  surface,  as 
defined  by  Equation  (5) 

Maximum  bubble  diameter  at  departure 
from  a  flat  surface 


ft 

ft 

ft 

ft 
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Volt 


e  Voltage  drop  across  resistance  R 

s  s 

during  calibration 

e  Voltage  drop  across  resistance  R 

during  calibration 

g  Earth  acceleration 

8  2 

g^  Conversion  factor,  4’18  x  10  ft/hr 

h„  Latent  heat  of  evaporation 

fg 


Volt 


ft/hr2 

lb  -ft/lb  -hr2 
m  f 

Btu/lb 

m 


I 

k 

K 

l 

l 

P 

i 

c 

q 


r 


r . 
i 


r 

o 


Current  flowing  through  the  heating  element 

Thermal  conductivity  of  liquid 

Thermal  conductivity  of  the  heating  element 

Thermal  conductivity  of  tube  at  the 

inner  wall  temperature 

Length  of  heating  element 

System  pressure 

Critical  pressure 

Heat  flux  at  the  surface  of  the 
heating  element 

Heat  generation  rate  in  the 
heating  element 

Radius 

Inner  radius  of  tube 

Outer  radius  of  the  heating  element 


amp 

Btu/hr-f t-F 
Btu/hr-f t-F 
Btu/hr-f t-F 

ft 

lbf/in2 
lb^/in2 
Btu/hr-f t2 

Btu/hr-f t3 

ft 

ft 

ft 


R  Electrical  resistance  ohm 

R.  Electrical  resistance  of  the  tube  at  ohm 

inner  wall  temperature 
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1 


m 


Electrical  resistance  of  the  tube  at 
mean  temperature 

Standard  electrical  resistance  used 
during  calibration 

Electrical  resistance  of  the  heating 
element  during  calibration 

Temperature  of  the  heating  element 
at  any  location 

Temperature  at  the  center  of  wire 


ohm 


ohm 


ohm 


T. 

1 


m 


w 

AT 

V 

a 


32 


v 

a 


Temperature  of  water  during  calibration 

Inner  tube  wall  temperature 

Mean  temperature  of  the  heating  element 

Saturation  temperature  of  water 

Outer  surface  temperature  of  the 
heating  element 

Wall  to  saturation  temperature  difference 

Voltage  drop  across  the  heating  element 

Temperature  coefficient  of  electrical 
resistivity 


F 

F 


Volt 

F-1 

.-1 


Temperature  coefficient  of  thermal  conductivity  F 
Resistivity  of  the  heating  element  at  32°F  ohm-ft 

3 


Density  of  liquid 
Density  of  vapor 
Surface  tension 


lb  /ft' 
m 


lb  /ft' 
m 


lb  /ft 
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Subscripts 

rt  During  the  room  temperature  calibration  test 

ht  During  boiling  at  atmospheric  pressure  calibration 
test 

expt  Experimental  value 
pred  Predicted  value 
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i 


ZFXX- 1 12-001 


AD-864  688 
TAB  70-T 


AD-876  732  abb 
Tab  7#- 1 
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